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Dual- Surface  Magnetic-Fieiu  Integral  Equation 

Solution  for  Bodies  of  Revolution 


1.  INTRODUCTION 


Conventional  solutions  of  the  magnetic-field  integral  equation  fail  to  produce  a  unique 
solution  at  frequencies  equd  to  the  resonant  frequencies  of  the  interior  cavity.*’*  These 
spurious  resonances  severely  corrupt  the  numerical  solution  of  the  magnetic-field  integral 
equation.  Yaghjian*  proved  that  the  original  magnetic-field  integral  equation  allows  spurious 
resonances  at  these  frequencies  because  it  does  not  restrict  the  tMgential  electric  field  to  zero 
on  the  surface  of  the  scatterer.  Tobin,  et  al.^  derived  the  dual-surface  magnetic-field  integral 

Received  for  Publication  1  May  1991 

*  Murray,  F.H,  (1931)  Conductors  in  an  electromagnetic  field,  Am.  J.  Math.,  63:275-288. 

^Maue,  A.W.  (1949)  On  the  formulation  of  a  general  scattering  problem  by  means  of  an  integral  equation, 
Zeitschrift  fur  Physik,  126(7/9)  ;601-618. 

^Yaghjian,  A.D.  (1981)  Augmented  electric  and  magnetic  field  integral  equation,  Radio  Science,  16:987- 
1001. 

^Tobin,  A.R.,  Yaghjian,  A.D.,  and  Bell,  M.M.  (1987)  Surface  integral  equations  for  multi-wavelength, 
arbitrarily  shaped,  perfectly  conducting  bodies,  Digest  of  the  National  Radio  Science  Meeting,  (URSI), 
Boulder,  CO. 
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equation  and  applied  it  to  a  three-dimensional,  multi-wavelength,  perfectly  conducting  body. 
This  numerical  solution  showed  that  the  dual-surface  magnetic-field  integral  equation  did 
indeed  eliminate  the  spurious  resonances  associated  with  the  interior  cavity  modes.  The 
solution  of  the  dual-surface  magnetic-field  integral  equation  first  appeared  in  Tobin,  et  al.*, 
and  the  proof  of  uniqueness  and  derivation  for  both  the  dual-surface  electric-field  integral 
equation  and  magnetic-field  integral  equation  will  appear®  in  March,  1991. 

Mautz  and  Harrington  derived  the  combined-field  integral  equation  solution®  to  elimi¬ 
nate  the  spurious  resonances  associated  with  the  interior  cavity  modes.  The  combined-field 
integral  equation,  as  the  name  implies,  combines  the  solution  of  the  magnetic-field  integral 
equation  with  the  electric-field  integral  equation.  The  combined-field  integral  equation’s 
main  drawback  is  that  it  involves  the  solution  of  both  the  magnetic-field  integral  equation 
and  the  electric-field  integral  equation,  therefore,  additional  programing  ability  and  computer 
run-time  are  required  compared  with  the  original  magnetic-field  integral  equation. 

One  application  of  the  combined-field  integral  equation  and  dual-surface  magnetic-field 
integral  equation  solutions  is  plane  wave  scattering  from  a  body  of  revolution.  As  an  alter¬ 
native  to  the  combined-field  integral  equation,  a  body  of  revolution  solution  was  formulated 
using  the  dual-surface  magnetic-field  integral  equation  to  determine  if  it  was  advantageous 
to  use  because  of  increased  accuracy  or  decreased  computer  run-time.  The  body  of  revo¬ 
lution  formulation  of  the  dual-surface  magnetic-field  integral  equation  solution  parallels  the 
body  of  revolution  formulation  in  Mautz  et  aJ.®  However,  the  numerical  solution  to  the  body 
of  revolution  dual-surface  magnetic-field  integral  equation  uses  a  single  pulse  and  impulse 
approximation  for  basis  and  testing  functions  instead  of  the  triangular  basis  and  testing  func¬ 
tions  used  in  Mautz  et  al.®  The  combined-field  integral  equation  solution  requires  a  more 
complicated  set  of  basis  and  testing  functions  due  to  the  inclusion  of  the  electric-field  inte¬ 
gral  equation.  The  different  basis  and  testing  fimctions  are  necessary  to  accurately  model 
the  derivative  of  the  siurface  current  in  the  electric-field  integral  equation.  This  derivative 
is  not  present  in  the  magnetic-field  integral  equation.  Also,  summations  for  the  integration 
use  a  simple  rectangular  rule  rather  than  the  more  complex  Gaussian  quadrature  integration 
scheme. 


^Woodworth,  M.B.,  and  Yagbjian,  A.D.  (1989)  Derivation,  application  and  conjugate  gradient  solution 
of  the  dual-surface  integral  equations  for  three-dimensional,  multi-wavelength  perfect  conductors,  PIERS-5: 
Applications  of  the  Conjugate  Gradient  Method  to  Electromagnetic  and  Signal  Analysis,  Sarkar,  T.K.  and 
Kong,  J.A.,  eds.,  ESsevier;  also  RADC  TR-89-142,  Hanscom  AFB,  MA,  ADA229076. 

^Mautz,  J.R.,  and  Harrington,  R.F.  (1987)  H-held,  and  oombined-held  solutions  for  bodies  of  revolution 
Arch.  Elecktron.  Ubertragungstedt.  (Electron.  Commun.),  32(4):157-164. 
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Results  from  the  dual-surface  magnetic-field  integral  equation  r  nd  combined-field  integral 
equation  body  of  revolution  solutions  are  compared  using  a  perfectly  conducting  sphere. 
The  accuracy  of  the  solutions  and  computer  run-times  are  also  compared  to  investigate  the 
benefits  and  pitfalls  of  both  solutions.  Far-field  results  for  the  truncated  and  flat-back  cones, 
as  well  as  the  truncated  circular  cylinder,  are  compared  with  the  results  of  a  wire-mesh 
computer  code  obtained  from  Ohio  State  University.^  The  predicted  radar  cross  section  of 
the  truncated  circular  e'dinder  is  also  compared  with  measured  data  over  a  wide  range  of 
aspect  angles. 


2.  STATEMENT  OF  PROBLEM 


It  is  desired  to  determine  the  scattered  far  field  from  a  perfectly  conducting  body  of 
revolution  excited  by  a  plane  wave.  Figure  1  shows  the  geometry  and  the  coordinate  system 
used  for  the  body  of  revolution,  where  p,  and  z  are  the  qylindrical  coordinates  and  t  and 
4>  form  an  orthogonal  curvilinear  coordinate  system  on  the  body  of  revolution  surface  S.  Ut 
and  are  orthogonal  unit  vectors  in  the  t  and  <f>  directions,  respectively.  Figure  2  shows 
the  coordinate  system  for  plane  wave  scattering.  The  propagation  vector  for  the  incident 
plane  wave  is  kt,  0t  defines  the  the  transmit  angle,  6r  and  <t>r  are  the  receiver^spherical  angles 
defining  the  direction  along  which  the  scattered  far  field  is  observed,  auid  fcr  is  the  recover 
propagation  vector,  which  points  from  the  receiver  location  to  the  origin.  Note  that  kt  is 
limited  to  the  xz-plane  {<f>t  =  0).  Figure  2  shows  that  uj,  wj,  and  are  unit  vectors  in 
the  6t,  y,  Or,  and  <f>r  directions  respectively.  The  notation  follows  basically  that  of  Mautz 
and  Harrington.® 


^Wang,  C.W.  (1983)  Scattering  From  Rotationally  Symmentric  Conducting  Body,  Ohio  State  University 
Electroscience  Laboratory  Tedi.  Rep.  714614-2,  Columbus,  Ohio. 
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Considered  scipaiately  is  an  incident  plane  wave  (//^)  defined  by 

(k,  X  ll'g)(  (1) 

and  an  incident  |)lane  \vavi‘  (//^)  flefined  by 

li,„r  (k,  X  u[,),  (2) 

where  Hinc  is  the  incident  inaffiietie  field,  f  is  the  position  vector  from  the  origin,  k  is  the  free 
space  wavenunibei .  and  //  is  the  free  s|)aee  intrinsic  impedance  (using  time  dependence), 
'rite;  incident  plane'  wave  gives  rise  to  electric  currents  directed  in  both  the  and  t  directions 
on  the  surface  S  as  well  as  scattered  far  fields  in  the  Or  and  (t>T  directions. 


3.  MAGNETIC-FIELD  INTEGRAL  EQUATION 


3.1  Formulation 


l  iguie  ilhistr.iies  an  aibitrarily  shai)ed  body  with  a  plane  wave  incident  on  the  body. 
I>\  dclinilioii  of  the  sc  atleicd  electric  and  inagnetic  fields, 

/•;(/■')  /v'.,,  t  i-\nr  (‘I) 

and 

ihr)  /l.sr  t  (4) 

I'oi  the  position  vector  r  not  on  the  .surface  5 

iLr(r)  -V  X  .-Cl/'),  (5) 

II.. 

where 

—  /  -4— (('•) 

4-  Is  r  -  r 


■\.Ar) 


X 


Figure  3.  Geometry  of  Perfect  Conductor  for  Plane  Wave  Scattering 
Substituting  for  in  Eq.  (5),  the  scattered  magnetic  field  becomes 


^u:(rl  =^^1^  J{f')  X  VVd5', 


where 

-}k\r-r\ 

For  r  inside  the  sphere,  E{f)  and  H(f)  equal  zero  from  the  boundary  conditions  at  the 
srirface,  and  thus 

HUn  ==  -ffincif).  (9) 

Substituting  from  Elq.  (9)  into  Eq.  (7)  we  obtain 

-//inc(r)  =  ^  f  J(f')  X  V'ifdS*  (f  inside  S).  (10) 

47r  Js 

As  r  approaches  S  from  inside  the  surface, 


/(f*)  X  VVd5'  =  £  Jif*)  X  VVdS'  +  ^  /(O  X  V'V»d5', 
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where  §  denotes  the  principal-value  surface  integral  evaluated  by  excluding  the  singular 
point  with  a  circular  "principal  area”  ^  and  f  is  the  integration  over  the  small  circular  patch 
as  r  approaches  5,  i.e.  as  e/6  — >  0  in  Figure  4.  At  r  near  r',  J(r)  «  J{r').  Therefore, 


^  J{f')  X  V'rl>dS'  =  /(r )  X  ^  Vxl>dS\ 

(12) 

Figure;  4  shows  that  for  the  small  patch  where  r  f' 

\r-f'\  =  p'^-\-e^  =  R'\ 

(Kl) 

and  therefore 

j  v'xi,ds'  =  j 

(14) 

where 

1 

II 

> 

(15) 

and 

R'  =  R'^p'  +  ITJ'  +  Kz'. 

(16) 

Substituting  for  e , 

R'  P'  n'j.  ^ 

-  R’" Ti’’ ■ 

(17) 

Evaluating  Eq.  (14)  with  Eq.  (17)  leads  to 

^  V  '■il)dS'  =  — 27rn, 

(18) 

and  therefore 

J{r)  X  ^  VVd5'  =  -J{f}  X  27rn. 

(19) 

Making  these  substitutions  in  Eq.  (10)  and  taking  n  cross 
field  integral  equation 

both  sides  we  obtain  the  magnetic- 

h  X  =  n  X  /  /  /(r')  x  V'^S' 

47r  ys  2 

(20) 
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Figure  4.  Geometry  of  Singularity  Integration 

This  derivation  was  for  a  circular  "principal  area”,  but  Yaghjian®  shows  that  this  deriva¬ 
tion  holds  for  any  shaped  “principal  area”  excluding  the  singularity  in  the  integrand  of  Eq. 
(20).  Using  the  vector  relationship 


/(O  X  w  =  -  y  jY), 


Eq.  (20)  becomes 
n  X  = 


8 


for  r  on  S.  The  outward  normal  n  is  defined  asn  =  u^xuf  Using  the  notation  in  Reference 
6  and  letting 

/(O  =  <},')  +  (23) 

where  u\  and  are  unit  vectors  in  the  if  and  <f>'  directions  respectively,  Eq.  (22)  expands  to 


Ut 


I  — ^  ^  lo  ~  P) 


—  {z'  —  z)  sin  xjf)  cos  4>'  —  2p  cos  vf  sin^(— )] 

2 

+  p'dt'  d(f>'GJ*{t',  <p'  +  <f>){z'  -  z)  sin  <f>^ 

+  I  +  ^P'^^  Jq  d<j/GJP{t\  <f>'  +  ^)(p'  sin  v  cos  v' 

— p  sin  t/  cos  u  —  (/  —  ^)  sin  v  sin  i/)  sin  0' 

j  <t>'de  jf"'  d<f>'GJ*ie,  4!  +  mkp'  -  p)  cos  u 

-(^-.)sin„)c<«^'  +  2,'cos«.„»(f))}=nx//,„, 


where 


and 


R=  yj{p-  p/y  +  (z  -  z'Y  +  4pp'  sin^Cy ). 


(24) 


(25) 


(26) 


In  Elq.  (24),  both  n  and  kt  are  to  be  evaluated  on  5  at  t  and  0.  The  angle  v  between  Ut  and 
the  z-axis  is  positive  when  tit  points  away  from  the  z-axis  and  is  negative  when  tit  points 
toward  the  z-axis.  The  primed  variables  are  functions  of  t'  and  therefore  are  evaluated  at  t'. 
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We  want  to  take  advantage  of  the  geometry  of  the  body  of  revolution  by  separating  the 
<f>  and  t  dependence  of  Jt  and  J^.  To  accomplish  this,  Jt  and  are  expanded  in  a  Fourier 
series  in  <f>: 


OO 


ns— OO 

(27) 

ns— OO 

(28) 

ns— OO 

(29) 

4/ +<!,)=  f;  iivy •<***>. 

(30) 

ns— oo 


The  Ut  and  comp)onent8  of  Eq.  (24)  will  be  evaluated  given  am  H*  polarized  incident 
wave.  For  the  Ut  component,  substituting  for  J‘(t,  <i>),  J*{t,  4),  (j/ +^),  and  4>' -^4) 

in  Eq.  (24)  yields 

f:  E 

ns— OO  ns— OO 

^K2{t\4')  £  =  (31) 

ns— OO 

where  Ki{t,tf,4')  and  K2{t,t',4')  are  the  known  portions  of  Eq.  (24)  that  depend  on  t,t', 
and  4'.  The  functions  of  t',  and  are  the  unknown  Fourier  coefficients.  The 

expression  on  the  right  side  of  Elq.  (31)  is  the  Ut  component  of  n  x  Hinc  for  the  H*  polarized 
incident  wave. 

Since  we  have  an  in  each  term  on  the  left  side  of  Eq.  (31),  we  will  multiply  both 
sides  by  and  integrate  both  sides  from  — tt  to  tt  with  respect  to  4-  The  integral  on  the 
left  side  of  Eq.  (31)  is  zero  if  m  ^  n  and  ir  if  m  =  n  (applying  orthogonality).  Evaluating 
the  right  side  of  Eq.  (24)  for  the  component  of  the  H*  polarized  incident  wave  where 

kff=—k{pcos4sindt  +  zcoadt),  (32) 

we  have 

r  [-fccoe^*<'’“'*""®'+*“'®‘^e--*'^]  d4.  (33) 

Jo 
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\\'(“  can  evaluate  this  inte'j;ial  usiiiff  liessi'l  hinelions  with 


/  j  r  r<  >s  0 

'i.:rL  ' 


for  m  =  0,±1.±2,. ..  (M) 

rh(-  ii'j,ht  side  of  1\(|.  (31)  using  Ecj.  (34)  ecjuals 

--kf  j"'  '\J„^^[kps\n(J^)  -  ./„_,(/i;/9sin6>t)].  (da) 

go  through  (he  same  substitution  as  for  the  ?/,  component  of  tlu;  77*^  polarized  incident 
wave  for  .7'(/.o).  ■l''{l.o).  +  (f)),  and  .7''’(<',</)'  f  0)  in  Ecj.  (24).  The  integral  on  the 

right  sidt'  of  luj.  (24)  is 


/'  |7sinosin  r( 

./(I 


I'sing  l'.'(|.  (.31)  redui  es  this  integral  to 

rA'Mii  (7 y''|./„_i(A'psin  (7,)  I  ./„  |  i(A:psinf7,)|.  (.37) 

When  we  u.se  th(!  saiiu;  approach  for  the  77®  polarized  incident  wave,  the  fp  component  is 
rA-eo.s77,r^'''""'®'y'*|.7„_i(A:psin6/e)  +  .J,n\{kps\nOt)\,  (38) 

and  the  ii^  conii)onent  is 

Ay }2rr  cos  r  kni  Oij  J,i{kp  Ot)  — 

sin  reostArr/'''' ‘(.7,,^  i(A-psiuf7<)  -  ./„+i  (Arpsin  fA,))].  (39) 

Ih'uiiiing  1-.<|.  (21).  we  find  that  the  tp  eomjManent  of  the  H’^  incident  wave  is 

.7„(U  ,  ^  f  I  1,1  I  1.1  111,:,’  I  !  \  I 


—  (;'  —  r)  sin  c')  cos()f)'  —  2pcos?7'.sin^(^)] 

'  y;  I  (l0'i  ‘"^''Ci{z  -  c)sin0' 


; a  '  ""  j.7„,  i(^'P^^iii(7,)  -  ./„_,(^7^sin(7,)l. 
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and  the  component  of  the  H*  incident  wave  is 

^  +  —p'dt'jnilf)  f  sin V cost/ 

2  47r  JQ 

— psint/coeu  —  (t/  —  z)8ini;8int/)sin^' 

j ^dHjtV)  £'  -  p)<x,sv 

±i 

— (/  —  ^)  sin  v)  cos  ^  +  2//  cos  t;  sin*(— )] 

=  ^sint;e’*'*“'®‘[;"“^J„+i(A:psin0t)  -/•J„_,(*:psin0|)]. 
For  the  /f*  incident  plane  wave,  the  tit  component  is 
^  ^  /  ^d£iUC)  d4.'^GW  -  p)  cos 


-(/  -  z)sint/)co8^  -  2poostf'sm*(^)l 

d4>^'^G{zf  -  z)  sin  (^' 

=  -^oos^,e''“'“®‘j"[yn+i(A:p8in^t)  +  yn-i(*;^sin^t)]. 
and  the  ti^  component, 


^  +  ^p'dt'j‘(0  /*'d0V"^G(0'sinvcost/ 

2  47r  ^0 

— p  sin  t/  cos  u  —  (/  —  z)  sin  u  sin  v')  sin 

J  ^  d^V’‘^G[((p'  -  p)  cos  v 

±9 

-{zf  -  z)  sin  v)  cos  ^'  +  2p'  cos  V  sin^(— )] 


(41) 


(42) 
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_  ^gjfcz oo« 9,  [2^  cos  ^  sin  OJ  Jn{kp  sin  6t) 


—  sin  V  cos  ( Jn+i  (fcp  sin  0^  —  J„+i  (Arp  sin  0t))] .  (43) 

The  t,  t'  and  <j/  integrations  can  be  converted  to  direct  summations.  Dividing  the  t  and  f 
integrals  into  discrete  segments,  we  evaluate  the  integrand  at  the  center  of  each  segment 
and  then  sum  over  the  number  of  t-  and  t'-segments.  The  <)>'  integral  is  divided  into  small 
patches,  again  approximating  the  integrand  at  the  center  of  each  patch.  In  moment  method 
terminology,  this  equates  to  using  pulse  basis  functions  and  unit  impulse  testing  functions. 
The  well-behaved  condition  of  the  magnetic-field  integral  equation  allows  such  simple  basis 
and  testing  fimctions  to  be  used. 

Rewriting  Eqs.  (40)  and  (41)  with  the  i!  and  integrations  represented  as  summations, 
we  have 

+  A1  f;  +  At  ^  =  Sl(t,A) 

m=l  m=l 

mjii  m^l 


for  /  =  1,2,...M,  (44) 


and 


^  +  At  f;  u)jL  +  At  x;  e.) 

m=l  m=l 

mfil  mjtl 


for  /  =  l,2,...Af,  (45) 

where  n  is  the  mode  number,  I  is  the  number  of  the  t-segment,  and  m  the  segment  number 
for  the  t/  integration.  The  first  superscript  indicates  the  direction  of  the  current  and  the 
second  superscript  indicates  the  vector  component.  M  is  the  maximum  number  of  t-  and 
t'-segments,  and  At  is  the  length  of  the  f-segment. 

For  a  given  t{l)  and  t'{m)  the  only  unknown  quantities  in  Eqs.  (44)  and  (45)  are  and 
j*.  Thus  for  a  given  n,  t,  and  t\  Eqs.  (44)  and  (45)  can  be  represented  as 
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■  5. 

. . . 

Ktt 

•••  ^Ihi 

’  j\  ■ 

■  s\  ■ 

k'u, 

•  •  • 

••• 

iL 

4 

Kt! 

^11 

...  A 

if 

St 

••  J 

.il. 

(46) 


where  K  and  S  are  known  functions  of  and  t  and  6t,  respectively  and  the  fs  are  the 

unknown  Fourier  series  coefficients  of  the  current. 

The  superscripts  of  the  j  and  S  variables  indicate  the  vector  component,  and  the  sub¬ 
scripts  indicate  the  segment.  This  square  matrix  equation  can  be  solved  for  and  for 
each  n,  and  the  currents  can  be  found  from  j*  and  j*  by  using  Eqs.  (27)  and  (^). 

Evaluating  the  <f/  integration  in  Eqs.  (40)  and  (43)  ai  t  =  t'  (p  =  f/  and  z  =  z')  yields 
R  =  2psin(<^/2).  Therefore,  a,t  t  =  tf,  the  phi  integration  does  not  converge.  This  lade 
of  convergence  is  not  due  to  the  magnetic-field  integral  equation  but  to  our  separation  of 
the  ^  and  tf  integrations.  To  obtain  convergence  at  t  =  f',  one  must  evaluate  the  <f>'lf  area 
integration  more  carefully.  In  Appendix  A,  it  is  shown  that  the  <//  integration  can  be  done 
properly,  when  the  (f/  integration  is  separated  from  the  If  integration,  if  the  modified  Green’s 
function 


replaces  G  in  Eq.  (25),  where 

H:  =  v'(A1/2)>  + 

(48) 

R.  =  v/(A1/2)V4  + 

(49) 

and  At  is  the  size  of  the  t-segment. 
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3.2  Far  Field 


To  obtain  the  scattered  magnetic  far  field,  let  r  ^  r'  in  Eq.  (7)  to  get 


Hscin  =  - 

47rr  Js 

(50) 

Noting  that 

f  X  u^  =  - 

■cos{<j>  —  (p')ug  +  cos^sin(<^  —  (i>')u4, 

(51) 

r  X  Ut  =  sin(<5!)  —  (p')  sin  v'ug  +  cos  6  cos{<f>  —  sin  v'u^ 

(52) 

~i> 

1 

II 

sin  Or  cos{<j)  —  (p!)  +  z' cos  B 

(53) 

Eq.  (50)  can  be  written  as 

HLirl  = 

Jo  /.  k‘(^'.<')(-C03(^ -,>')) 

+ 

J‘(^',t')(sin(^  -  <^')sini;')] 

(54) 

and 

Htli^  = 

47rr  Jo  Jo  (p)cosB) 

+ 

J*{(p',  <')(  —  sin^cos  u'  +  cos{<p  —  ^')  cos  0  sin  u')j 

^jk{p'  sin  e  co»{<t)-<t>')+i'  cos 

(55) 

where,  according  to  Eqs.  (27)  and  (28) 

w,o  = 

oo 

n=— oo 

(56) 

1 

8 

;s 

(57) 

15 


16 


For  the  plane  wave  incident  along  the  positive  z-axis,  only  the  ±  1  modes  are  needed 
in  the  Fourier  series  to  represent  and  For  an  incident  plane  wave  off  the  z-axis,  the 
number  of  modes  n  was  chosen  as 

n=  (l.lA:rsin0t  +  2)  (63) 

where  r  is  the  maximum  radius  of  the  body  of  revolution  measured  from  the  z-axis.®  The 
reactive  fields  of  the  scatterer  are  extremely  small  for  n  >  kr.  Assuming  the  number  of  far- 
field  fluctuations  per  solid  angle  are  roughly  the  same  in  all  directions,  the  far- field  variations 
with  (j)  will  be  less  near  the  z-axis.  Therefore,  the  number  of  modes  needed  to  expand  the  field 
will  become  less  as  one  approaches  the  z-axis,  since  the  higher  order  modes  which  represent 
the  more  rapidly  varying  parts  of  the  field,  will  not  be  needed  to  expand  the  slower  variations 
with  4>.  This  demonstrates  that  the  maximum  number  of  modes  needed  to  represent  the  !  ax 
field  at  the  angle  6t  is  approximately  ka  sin  6t.  The  addition  of  the  extra  two  modes  was 
determined  by  the  examination  of  results  of  the  code  for  different  number  of  modes. 

Typically  five  to  ten  segments  per  wavelength  are  used  for  both  the  number  of  t-segments 
and  phi-segments.  The  criterion  for  the  number  of  phi-segments  is  based  on  having  square 
patches  for  the  phi  integration.  Other  shaped  patches  were  used,  but  the  results  proved  most 
accurate  for  square  or  approximately  square  patches. 


3.3  Numerical  Results 


Figure  5  is  a  plot  of  the  normalized  radar  cross  section  a/Tro^  in  the  backseat tering 
direction  from  a  sphere  computed  using  the  magnetic-field  integral  equation  and  the  Mie 
scries^’ with  the  incident  wave  along  the  positive  z-axis.  The  Mie  scries  is  an  exact  an¬ 
alytical  solution  for  the  scattering  from  spheres  only  and  cannot  be  used  to  validate  the 
dual-surface  magnetic-field  integral  equation  solution  for  other  geometries. 


*Yaghjian,  A.D.  (1977)  Near  Field  Antenna  Measurements  on  A  Cylindrical  Surface:  A  Source  Scattering- 
matrix  Formulation,  NBS  Technical  Note  696. 

^DiBeneditto,  J.P.  (1984)Bistatic  Scattering  from  Conducting  Calibration  Spheres,  RADC-TR-84-93, 
Hanscom  AFB,  MA,  ADA154173. 

*°Stratton,  J.A.  (1941)  Electromagnetic  Theory,  McGraw-Hill,  New  York. 


17 


Figure  5  shows  that  the  magnetic-field  integral  equation  solution  contains  spurious  reso¬ 
nances  at  the  interior  cavity  resonance  frequencies  of  2.744,  3.870,  4.973,  5.763,  6.062,  6.988, 
7.140,  8.183,  8.211,  9.275,  and  9.356."  The  limitations  of  the  magnetic-field  integral  equa¬ 
tion  are  clearly  demonstrated  by  Figure  5. 


4.  DUAL-SURFACE  MAGNETIC-FIELD  INTEGRAL  EQUATION 


4.1  Formulation 


The  following  derivation  and  proof  of  uniqueness  of  the  dual-surface  magnetic-field  in¬ 
tegral  equation  closely  follows  the  work  on  the  dual-surface  magnetic-field  integral  equation 
presented  by  Yaghjian.®  The  dual-surface  magnetic-field  integial  equation  can  be  derived 
by  beginning  with  Eq.  (10),  the  “interior”  or  “extended”  magnetic-field  integral  equation, 
which  is  rewritten  here 

-Hincin  =  f  Jif')  X  VipdS'  (f  inside  5).  (64) 

47r  Js 

The  cuTTent  J{f)  in  Eq.  (64)  is  uniquely  determined  at  every  frequency  if  Eq.  (64)  is  satisfied 
for  all  r  inside  5.*^ 


^‘Harrington,  R.F.  (1961)  Time-Hanruxtic  Electromagnetic  Fields,  McGraw-Hill,  New  York. 
‘^Waterman,  P.C.  (1965)  Matrix  formulation  of  electromagnetic  scattering,  Proc.  IEEE,  53:805-812. 
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DIL 


ka 


Figure  5.  Backscattering  Cross  Section  Versus  ka  of  a  Perfectly  Conducting  Sphere  Com¬ 
puted  with  the  Original  Magnetic-Field  Integral  Equation  and  the  Mie  Series 
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Figure  6  shows  the  surface  S  and  a  surface  Sg  pairallel  to,  and  some  distance  6  >  0, 
inside  of  S. 

By  adding  an  cross  Exj.  (64)  at  points  f  on  Sg  to  the  corresponding  points  on  S  in  the 
original  magnetic-field  integral  equation,  one  obtains  the  dual-surface  magnetic-field  integral 
equation:® 


n  X  Ho{r^  =  -  n  x 


VMry)dS‘ 


where 


Ho^  =  Hincif)  +  Sn) 

ipoify)  =  ipify)  -f  aip{f-  Shy) 
and,  as  usual 


(65) 


(66) 

(67) 


(68) 


Figure  6.  Geometry  of  Perfect  Conductor  for  Dual-Surface  Magnetic-Field  Integral  Equation 
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The  dual-surface  magnetic-field  integral  equation  will  have  a  unique  solution  for  J  at  all 
real  frequencies  if  a  is  imaginary  and  the  positive  real  constant  S  is  less  than  approximately 
A/2. 

By  considering  the  fields  radiated  by  the  solution  currents,  uniqueness  of  solution  can  be 
shown  for  the  dual-surface  magnetic-field  integral  equation.  ^s(f)  will  be  the  magnetic  field 
radiated  by  the  solution  J 

H,{f)  =  ^  /  /(r^)  X  V'rl)dS'  (f  not  in  S).  (69) 

47r  Js 

If  J  is  the  unique  scattering  current  for  Eq.  (64),  then  is  the  correct  scattered  field 

for  all  r  not  in  the  surface  current.  Because  it  is  not  known  that  J  is  the  correct  unique 
solution,  Eq.  (64)  only  defines  an  unknown  magnetic  field 

Taking  the  curl  of  the  curl  of  Eq.  (64)  shows  that  this  unknown  magnetic  field  satisfies 
the  homogeneous  vector  wave  equation  for  all  r  not  in  the  surface  current  J,  that  is, 

V  xV  X  H,  -  k^Ht  =  0  (r  not  in  5).  (70) 

Let  f  approach  S  in  Eq.  (64)  from  the  inside  of  S  and  using  the  principal  value  formula^ 

/  Vtl>dS'  =  /  VxPdS'  -  27rn,  (71) 

Js(r^S)  Js 

as  explained  above  in  the  derivation  of  the  magnetic-field  integred  equation,  Eq.  (64)  can  be 
written  as 

=  ^l  J(r  0  X  VV(r,f')d5'  -  fi  x  ^  (f  on  5),  (72) 

where  f_  indicates  the  field  evaluated  just  inside  the  surface  current.  Since  Eq.  (64)  is  true 
for  all  r  inside  5,  Ht  on  the  parallel  surface  Ss  can  be  written  as 

H, =  ^  ^V(r  -  6h)dS'  (f  on  5) .  (73) 

Multiply  Eq.  (73)  by  a,  add  to  Elq.  (72),  and  then  take  h  cross  the  results  to  obtain 

n  X  [/7.(f_)  -I-  aH,{f-  6n)]  =  n  x  J{f')  x  V'tito(r,r')d5'  -  ^ 

(f  on  5).  (74) 
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(75) 


Comparing  Eqs.  (74)  and  (65),  which  J  must  satisfy,  shows 

n  X  [Hg{f-)  +  Hincif)  +  aH,{f-  6h)  +  aHindf-  5n)]  =  0  (f  on  S). 
The  incident  magnetic  field  must  also  satisfy  the  vector  wave  equation 


V  X  V  X  Mine  ~  k^Hinc  =  0.  (76) 

Adding  Eq.  (70)  to  Elq  (76)  and  rewriting  Eq.  (76)  to  arrive  at  the  interior  boundary  value 
problem 

VxVxH-k‘^H  =  Q  (r  inside  S),  (77) 

n  X  [^(r)  +  aH(f—  6n)]  =  0  (f  — >  5  from  inside),  (78) 

where  H(r)  is  the  total  magnetic  field  given  by 

(79) 


Now  it  needs  to  be  shown  that  the  boundary  value  problem  of  Eqs.  (77)  and  (78)  has 
only  the  trivial  solution,  H{i^  =  0,  for  the  total  field  throughout  the  volume  enclosed  by 
the  surface  5,  provided  the  constant  a  is  imaginary  and  the  positive  real  constant  6  is 
smaller  than  approximately  A/2.  Rewrite  the  boundary  condition  Eq.  (78)  explicitly  for  the 


magnetic  field  tangent  to  the  the  surface  S  to  get, 

Htif)  +  o^t(f  —  6n)  =  0  {f  S  from  inside).  (80) 

Rewrite  the  tangential  magnetic  fields  in  Eq.  (80)  in  terms  of  a  magnitude  and  phase 

fitif)  =  |^i(f)|e'^  if-*  S  from  inside),  (81) 

Ht{r-  6h)  =  [|M(r)(  +  A/f,(r, d)]  (f-  5  from  inside),  (82) 

where  A//t  and  are  the  differences  between  the  magnitudes  and  the  phases  of  Ht{f)  and 
6h).  Now  using  Eqs.  (81)  and  (82)  in  Eq.  (80)  results  in 

\Ht\  -i-  a{\Ht\  +  AHt) (cos  A^  +  y  sin  A^)  =  0.  (83) 

Since  \Ht\,  A^t,  and  A^  are  real  numbers,  by  letting  the  constant  a  be  an  imaginary  number 
(Joj),  we  can  equate  the  real  and  imaginary  parts  of  Eq.  (83)  separately  to  give 

\Ht\-ai(\Ht\  +  AHt)smA(l>  =  0,  (84) 
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and 


+  A^t)cos  A<;6  =  0.  (85) 

For  small  6,  A^  will  be  small  (much  less  than  ±90°),  and  thus  Eqs.  (84)  and  (85)  imply 
that  |/7t|  =  0  and  A  fit  =  0,  or 

fftii^  =  0  (f  S  from  inside)  (86) 

^f(r  — 6n)  =  0  (r —»  S  from  inside).  (87) 

This  says  that  when  the  constant  a  is  chosen  imaginary  and  6  is  not  large,  the  two  separate 
tangential  fields  in  the  boundary  condition  are  zero,  meaning  that  the  tangential  magnetic 
field  on  both  surfaces  is  zero. 

The  boimdary  condition  of  Eq.  (86)  restricts  the  nonzero  solutions  of  Eq.  (77)  to  the 
resonant  modes  of  the  cavity  formed  by  a  perfectly  magnetically  conducting  surface  S.  These 
modes  form  standing  waves  within  the  cavity  with  magnetic  and  electric  fields  that  can  be 
chosen  real  and  imaginary,  respectively.*^  The  tangential  magnetic  field  near  the  surface  can 
be  expressed  approximately  as 

Tn)  «  A(ft,  Tn)  sin  7r„,  (88) 

where  (r*,r„)  are  the  coordinates  tangent  and  normal  to  the  surface  5,  7  is  a  pKJsitive  real 
propagation  constant  with  a  value  equal  to  or  less  than  the  free  space  constant  k  ,  and  the 
amplitude  >l(f„r„)  varies  with  r„  slowly  compared  to  the  variation  of  sin7r„.  If  r„  equals 
zero  on  5,  the  boxmdary  condition  of  Eq.  (86)  is  satisfied  by  Elq.  (88).  The  boundary 
condition  of  Eq.  (87)  applied  to  Eq.  (88)  requires  that 

76  =  mir  (89) 

for  m  equal  to  a  positive  integer.  (It  is  assumed  is  that  there  will  be  some  portion  of  S  where 
A  is  nonzero.)  Since  the  maximum  value  of  7  is  A:  =  2ir/\,  Eq.  (89)  cannot  be  satisfied  for 

0<6<A/2.  (90) 

An  approximate  sign  is  included  in  Eq.  (90)  because  Eq.  (88)  is  an  approximate  expression 
for  the  standing  wave  field  near  the  surface. 

'^Borgnis,  F.E.  and  P^pas,  C.H.  (1958)  Electromagpietic  waveguides  and  resonators,  Encyclopedia  of 
Physics,  S.  Flugge,  ed.,  16:  Springer- Verlag,  Berlin. 
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To  summarize,  the  only  solution  to  Elqs.  (77)  and  (78)  for  a  imaginary  and  0  <  6  ;$  A/2 
is  the  trivial  solution,  l^(f)  =  Hincif)  +  =  0  throughout  the  volume  enclosed  by  S. 

Since  =  —V  x  H/jueo,  the  electric  field  E(f)  =  Emcif)  +  E,{f)  within  this  volume  is 
also  identically  zero.  It  can  be  proved^^  that  the  current  that  produces  the  negative  of  the 
incident  electromagnetic  fields  throughout  the  volume  enclosed  by  S  is  the  correct  rmique 
current  for  the  exterior  scattering  problem.  Since  this  unique  solution  is  derived  from  the 
solution  current  of  Elq.  (65),  the  dual-surface  magnetic-field  integral  equation  has  a  unique 
solution. 

The  existing  magnetic-field  integral  equation  solution  is  easily  modified  to  obtain  the 
dual-surface  magnetic-field  integral  equation.  The  values  of  Hinc  and  V'  are  now  calculated 
at  both  the  outer  surface  and  an  irmer  surface  some  distance  6  inside.  This  requires  that  the 
computer  program  be  modified  to  simply  evaluate  these  functions  at  r  —  6n  as  well  as  at  f. 

For  the  body  of  revolution  solution,  a  was  chosen  to  be  j  (\/-T)  to  weight  the  field  on 
S  and  Ss  equally  in  the  boundary  conditions  of  Eqs.  (77)  and  (78).  To  keep  the  surface  St 
about  an  equal  distance  between  the  two  critical  values,  6  —0  and  A/2,  6  was  selected  to 
equal  A/4.® 


4.2  Numerical  Results 


Figure  7  plots  the  normalized  radar  cross  section  o/tto*  in  the  backscattering  direction 
from  the  same  spherical  scatterer,  as  we  did  above  with  the  magnetic-field  integral  equation, 
using  the  dual-surface  magnetic-field  integral  equation  and  the  Mie  series  with  the  incident 
wave  along  the  positive  z-axis.  The  spurious  resonances  are  eliminated  by  use  of  the  dual- 
surface  magnetic-field  integral  equation,  and  there  is  excellent  agreement  out  to  ka  =  10. 
At  ka  =  10,  there  are  five  t-segments  per  wavelength,  and  as  this  number  of  t-segments  falls 
below  five  per  wavelength,  the  dual-surface  magnetic-field  integral  equation  solution  rapidly 
degrades.  This  is  expected  because  a  general  rule  of  thumb  for  moment  method  solutions 
for  the  magnetic-field  integral  equations  is  to  have  no  less  than  about  five  segments  per 
wavelength.® 
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Figure  7.  Backscattering  Cross  Section  Versus  ka  of  a  Perfectly  Conducting  Sphere  Com¬ 
puted  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  Mie  Series 
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Figure  8  shows  the  comparison  of  the  Mie  series,  dual-surface  magnetic-field  integral 
equation,  and  the  combined-field  integral  equation*  for  the  same  spherical  model  with  the 
incident  wave  along  the  p>ositive  z-axis.  The  dual-surface  magnetic-field  integral  equation 
and  combined-field  integral  equation  are  in  good  agreement  up  to  iba  =  5  (ten  segments 
per  wavelength)  where  the  combined-field  integral  equation  starts  to  deviate  from  the  Mie 
series  solution.  As  stated  previously,  the  dual-surface  magnetic-field  integral  equation  does 
not  show  a  significant  deviation  until  ka  =  10.  The  computer  run  times  for  the  dual-surface 
magnetic-field  integral  equation  and  combined-field  integral  equation  were  approximately 
equal  for  the  spherical  model  with  the  same  number  of  patches  per  wavelength.  However, 
the  dual-surface  magnetic-field  integral  equation  requires  less  segments  per  wavelength  for 
the  same  accuracy. 

The  dual-surface  magnetic-field  integral  equation  solution  was  computed  as  a  function  of 
bistatic  angle  for  a  sphere  of  iba  =  20  with  the  incident  field  along  the  positive  z-axis.  This 
result  was  compared  with  the  Mie  series  solution,  and  again  the  agreement  is  excellent  as 
seen  in  Figures  9  and  10. 

To  ensure  that  an  adequate  number  of  modes  are  being  calculated,  according  to  the 
criterion  of  Eq.  (63),  the  same  spherical  model  was  used  and  the  incident  field  was  moved  to 
the  positive  i-axis.  Figures  11  and  12  show  the  comparison  of  the  dual-surface  magnetic- field 
integral  equation  and  the  Mie  series.  Again,  the  agreement  is  excellent. 

To  verify  the  accuracy  of  the  solution  for  bodies  large  with  respect  to  wavelength,  the 
cross  section  of  a  sphere  with  ka  =  70  with  the  incident  field  coming  in  along  the  positive 
z-axis  was  also  computed.  These  results  are  summarized  in  Figrires  13  and  14.  Due  to  a 
limitation  on  computer  time,  however  the  case  of  the  same  sphere  with  the  incident  field 
along  the  positive  x-axis  was  not  computed.  An  attempt  was  also  made  to  obtain  a  result 
from  the  combined-field  integral  equation  with  the  incident  field  along  the  positive  x-axis 
for  a  sphere  with  ka  =  20,  but  this  was  unfruitful.  The  results  obtained  were  meaningless, 
but  the  reason  for  this  is  unclear. 

Although  the  sphere  results  show  that  the  dual-surface  magnetic-field  integral  equation 
solution  appears  to  work  well  for  a  sphere,  results  from  other  geometries  would  increase 
the  confidence  in  its  answers.  Wang^  plotted  the  results  of  plane  wave  scattering  from 
rotationally  symmetric  conducting  bodies  computed  from  a  wire-mesh  code.  A  copy  of 
the  code  was  obtained  and  several  cases  presented  in  Reference  7  were  run.  Figures  15-24 
compare  the  results  of  the  dual-surface  magnetic-field  integral  equation  and  the  Ohio  State 
code  applied  to  truncated  and  flat-back  cones  and  the  circular  cylinder.  The  radar  cross 
section  was  calculated  as  a  function  of  bistatic  smgle  for  incident  plane  waves  along  the 
positive  and  negative  z-axes,  respectively. 
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Figure  8.  Backscattering  Cross  Section  Versus  ka  of  a  Perfectly  Conducting  Sphere  Com¬ 
puted  with  the  Dual-Surface  Magnetic-Field  Integral  Elquation,  Combined-Field  Integral 
Equation,  and  the  Mie  Series 
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Figure  10.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Sphere 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  Mie  Series;  ka  = 
20,  dt  =  0,  H-Plane 
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Figure  11.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Sphere 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  Mie  Series;  ka  = 
20,  Ot  =  90,  E-Plane 
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Figure  12.  Raxiar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Sphere 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Elquation  suid  the  Mie  Series;  ka  = 
20,  $t  =  90,  H-Plane 
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Figure  13.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Sphere 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  Mie  Series;  ka  = 
70,  6t  =  0,  Emplane 
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Figure  14.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Sphere 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  Mie  Series;  ka  = 
70,  dt  =  0,  H-Plane 


Figures  15  and  16  show  no  signihcant  difference  in  the  results  of  the  dual-surface 
magnetic-field  integral  equation  and  OSU  code  for  the  circular  cylinder.  As  the  circulau’ 
cylinder’s  radar  cross  section  is  identical  for  the  incident  wave  along  the  positive  z-axis  and 
the  negative  z-axis  as  a  function  of  bistatic  angle,  only  the  case  of  the  incident  wave  along 
the  positive  z-axis  is  shown. 

Figures  17-20  show  the  results  of  the  truncated  flat-back  cone.  I  he  S-plane  results  for 
the  two  codes  again  show  no  significant  difference.  The  H-plane  results  have  slight  deviations 
of  up  to  1  dB  but  are  almost  identical  for  the  whole  range  of  bistatic  angles. 

The  results  for  the  flat-back  cone  are  shown  in  Figures  21-24.  The  maximum  difference 
in  the  results  for  the  E-plane  is  approximately  3dB.  This  occurs  for  the  incident  wave  along 
the  positive  z-axis  from  about  135  to  140  degrees  bistatic  angle.  The  E-plane  results  for  the 
incident  field  along  the  negative  z-axis  are  in  excellent  agreement  showing  differences  of  a 
fraction  of  a  dB.  The  difference  in  the  H-plane  results  for  the  flat-back  cone  are  the  opposite 
of  the  El-plane  results.  There  is  less  than  a  0.5  dB  difference  for  the  incident  field  along  the 
positive  z-axis  and  a  maximum  of  3  dB  difference  from  40  to  110  degrees  bistatic  angle  for 
the  incident  field  along  the  negative  z-axis. 

OveraU  there  is  good  agreement  between  the  two  codes.  Figure  25  shows  the  results  from 
the  OSU  code  for  a  circular  cylinder  for  backscatter  versus  ka  where  o  is  the  radius  and 
one-half  the  length  of  the  cylinder.  As  ka  increases,  both  the  length  and  diameter  of  the 
cylinder  change  appropriately.  It  is  seen  that  the  OSU  code  has  spurious  resonances  like  the 
magnetic-field  integral  equation  solution.  These  spurious  resonances  may  be  the  reason  for 
the  difference  in  the  results  of  the  two  codes. 

Another  reason  may  be  the  difference  method  used  to  approximate  the  currents  on  the 
body  with  the  OSU  code  using  wire-segments  and  the  dual-surface  magnetic-field  integral 
equation  solution  using  patches.  The  results  for  the  dual-surface  magnetic-field  integral 
equation  were  obtained  using  ten  segments  per  wavelength  and  the  OSU  code  with  11.8 
segments  per  wavelength.  Figure  26  shows  the  comparison  of  the  results  for  the  OSU  code 
for  the  flat-backed  cone  using  10.0  and  11.8  segments  per  wavelength.  There  is  a  1  dB 
difference  in  the  two  results  at  backscatter.  This  seems  to  indicate  that  the  dual-surface 
magnetic-field  integral  equation  has  increased  accuracy  over  the  OSU  code  for  the  same 
number  of  segments  per  wavelength  as  the  results  from  the  two  codes  show  little  difference 
if  the  OSU  code  uses  11.8  segments  per  wavelength. 

Figures  27  and  28  show  the  El-plane  and  H-plane  monostatic  scattering  from  a  circular 
cylinder  respectively.  The  measurements  were  performed  at  the  Rome  Laboratory’s  electro¬ 
magnetic  experimental  scattering  facility.  The  results  for  both  polarizations  show  a  difference 
of  approximately  1.5  dB  at  the  specular  angles.  Modifying  the  measured  values  by  the  dif¬ 
ference  of  the  physical  optics  approximation  and  the  measurements  at  the  specular  angles 
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would  produce  excellent  agreement  for  the  El-plajie  results.  A  larger  difference  in  the  H-plane 
results  would  be  due  to  experimental  error,  as  is  demonstrated  by  the  lack  of  symmetry  in 
the  measured  results.  The  experimental  error  for  the  measurements  is  ±  1.5  dB  and  the 
difference  between  the  measurements  and  predicted  values  are  within  these  bounds.  For 
monostatic  results,  the  body  of  revolution  code  must  compute  a  new  current  solution  at 
every  observation  angle.  The  results  for  the  circular  cylinder  with  a  radius  of  1.2675  A  and 
a  length  of  6  A  required  72  hours  of  CPU  time  on  a  VAX6000.  This  enormous  amount  of 
computer  time  prevented  the  comparison  of  the  dual-surface  magnetic-field  integral  equation 
with  other  measured  monostatic  results. 
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Figure  15.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Circular 
Cylinder  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  OSU 
Code;  6t  =  0,  E-Plane 
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Figure  16.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Circular 
Cylinder  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Elquation  and  the  OSU 
Code;  Bt  =  0,  H-Plane  „„ 


Figure  17.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Frustum 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  OSU  Code; 
6t  =  0,  EJ-Plane 


Figure  18.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Frustum 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Elquation  and  the  OSU  Code; 
6t  =  0,  H-Plane 
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Figure  19.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Frustum 
Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  OSU  Code; 
=  180,  E-Plane 
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Figure  20.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Prustum 
Computed  with  the  Dual-Surface  Magnetic-Field  IntegrEJ  Elquation  and  the  OSU  Code; 
6t—  180,  H-Plane 
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Figure  21.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Flat-Backed 
Cone  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  OSU  Code; 
dt  =  0,  EJ-Plane 
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Figure  22.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Flat-Backed 
Cone  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  OSU  Code; 
6t  =  0,  H-Plane 
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Figure  23.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Flat-Backed 
Cone  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  OSU  Code; 
Ot  =  180,  &Plane 
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Figure  24.  Radar  Cross  Section  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Flat-Backed 
Cone  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Equation  and  the  OSU  Code; 
6t  —  180,  H-Plane 


45 


(J/7Ta 


(j/(47rX^  )  (dB) 


Figure  26-  Radar  Cross  Scclion  Versus  Bistatic  Angle  of  a  Perfectly  Conducting  Flat-Backed 
Cone  Computed  with  the  Dual-Surfau^e  Magnetic-Field  Integral  Equation  (10.0  Segments  Per 
Wavelength)  and  the  OSU  Code  (10.0  and  ll.S^cgments  Per  Wavelength);  Ot  =  180,  E-Plane 
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Figure  27.  Backscattering  Cross  Section  Versus  Aspect  Angle  of  a  Perfectly  Conducting 
Circular  Cylinder  as  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Ekjuation 
and  Measurements;  E-Plane 


Figure  28.  Backscattering  Cross  Section  Versus  Aspect  Angle  of  a  Perfectly  Conducting 
Circular  Cylinder  as  Computed  with  the  Dual-Surface  Magnetic-Field  Integral  Elquation 
and  Measurements;  H-Plane 
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5.  CONCLUSIONS 


It  has  been  demonstrated  that  the  dual-surface  magnetic-field  integral  equation  eliminates 
the  spurious  resonances  associated  with  the  magnetic-field  integral  equation  for  plane  wave 
scattering  from  bodies  of  revolution.  It  was  an  extremely  simple  modification  to  convert  the 
original  magnetic-field  integral  equation  computer  program  to  the  dual-surface  magnetic- 
field  integral  equation  program.  As  the  magnetic-field  integral  equation  is  very  well  behaved, 
simple  pulse  basis  and  impulse  testing  functions  provide  excellent  results  for  as  few  as  five 
segments  per  wavelength. 

The  accuracy  of  the  dual-surface  magnetic-field  integral  equation  was  shown  to  be  equal 
or  better  than  the  accuracy  of  the  combined-field  integral  equation  for  a  perfectly  conducting 
sphere.  The  combined-field  integral  equation  incorporates  both  the  magnetic-field  integral 
equation  and  the  electric-field  integral  equation,  which  requires  more  complicated  basis  and 
testing  functions  for  accurate  results.  These  more  complicated  functions,  along  with  the 
necessity  to  incorporate  the  the  electric-field  integral  equation,  malce  the  combined-field 
integral  equation  solution  more  difficult  to  program  than  the  dual-surface  magnetic-field 
integral  equation  solution.  Even  with  the  more  complicated  basis  and  testing  functions, 
the  combined-field  integral  equation  started  to  decrease  significantly  in  accuracy  below  ten 
segments  per  wavelength.  This  was  not  true  for  the  dual-surface  magnetic-field  integral 
equation.  The  computer  run-times  of  the  two  solutions  were  approximately  equal  for  the 
perfectly  conducting  sphere  for  equal  number  of  segments  per  wavelength. 

The  dual-surface  magnetic-field  integral  equation  results  were  also  in  good  agreement 
with  results  produced  by  a  wire-mesh  body  of  revolution  code.^  The  results  were  identical 
for  a  circular  cylinder  but  deviated  by  as  much  as  3  dB  for  the  truncated  and  flat-back  cones. 
The  differences  in  the  two  results  may  be  due  to  the  different  methods  used  to  approximate 
the  current  on  a  solid  surface  and  the  fact  that  the  OSU  code  does  not  eliminate  the  spurious 
resonances. 

The  comparison  of  the  measurements  with  the  dual-surface  magnetic-field  integral  equa¬ 
tion  results  for  a  circular  cylinder  also  were  quite  good,  demonstrating  that  the  dual-surface 
magnetic-field  integral  equation  computes  accurate  results  over  a  wide  range  of  incident  an¬ 
gles.  The  difference  in  the  results  from  the  two  different  polarizations  were  within  the  limits 
of  experimental  error.  Due  to  the  excessive  computer  run  time  required  to  obtain  these 
results,  however,  no  other  measured  monostatic  results  were  compared  with  the  dual-surface 
magnetic-field  integral  equation  solution. 
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The  next  obvious  step  is  to  attempt  to  decrease  the  computer  run  time  of  the  dual-surface 
magnetic-field  integral  equation  solution  when  the  angle  of  incidence  is  far  from  the  2-axis  of 
revolution.  This  can  be  accomplished  by  using  a  fast  Fourier  transform  (FFT)  to  compute 
the  phi  integration  for  each  mode.  The  FFT  will  save  computer  run  time,  for  cases  with 
the  incident  field  off  the  2-axis,  of  about  21og2iV/JV  times  the  current  run  time  where  N  is 
approximately  kr  (see  Elq.  (63))  for  the  body  of  revolution  program.  The  drawback  to  the 
FFT  is  the  additional  memory  requirement.  This  could  be  overcome  by  using  direct  access 
files  on  disk  as  is  discussed  in  Woodworth.*^  By  using  direct  access  files,  however,  there  is 
the  added  difficulty  in  computer  programming,  a  greater  CPU  time  requirement,  and  the 
possibility  of  a  large  increase  in  input /output  time. 

The  addition  of  a  geometiy  package  allowing  the  user  to  input  a  variety  of  geometry  types 
would  greatly  aid  the  use  of  the  program.  At  the  present  time,  the  program  accepts  spheres, 
conespheres,  cylinders,  and  cylindrical  shaped  fiiistums  as  input.  The  code  was  originally 
written  to  accept  only  a  sphere  as  an  input  but  the  geometry  input  package  of  Reference  7 
was  modified  to  allow  the  other  shapes  to  be  used  as  input  to  the  code.  Other  shapes  would 
provide  additional  validation  of  the  dual-surface  magnetic-field  integral  equation  solution  as 
well  as  providing  useful  results  to  the  users  of  the  body  of  revolution  program. 


‘^Woodworth,  M.B.  (1988)  Large  Matrix  Solution  Techniques  Applied  to  an  Electromagnetic  Scattering 
Problem,  RADC-TR-8i268,Hanscom  AFB,  MA,  ADA206917. 
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Appendix  A 

Integration  Correction 


The  purpose  of  Appendix  A  is  to  obtain  Ekj.  (47).  For  this  integration  let  t'  —  t  =  I'  by 
letting  t  =  0.  As  shown  in  Figure  Al,  the  integration  of  the  Green’s  function  for  the  ring 
containing  the  self  patch  involves  the  integral 

(AO 

Jo 


113/2 


This  can  be  rewritten  as 

I  ^  7^/  072 _ , , 

Jo  ^  ^  JAt/2  {{p  ~  pfY  +  {z  -  z'y  +  4pp'sin^((;^/2)]'V2 


(A2) 


by  substituting  for  H  from  Eq.  (2(i).  For  a  fixed  <p'  the  t'  integiation  in  E<|.  (A2)  can  be 
written  ;is 
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At/2 

At/2 


1 


dt'. 


(A3) 
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For  </»'  small 


where 


Rc  =  p4>' 


and 


('  =  v/(p  -  +  (z- 

Therefore 


R  «  \/pV^  + 


and  so 


/,  =2 


I 


At/2 


1 


:dt'. 


For  a  small  fixed  4/  and  integrating  Ekj.  (A8),  one  obtains 


At/2 


P^4>'"  [V^(At/2)2+p20t2j 

Im)!  a  small  (ixcd  0'  and  At  a])proaching  zero, 


h  = 


At. 


The  same  result  will  he  obtained  for  a  small  fixed  <p'  from  the  following  etjuation 
At 


I 


8/r’  sin^  4)'/2 


(A4) 

(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(AlO) 

(All) 


h() 


This  \/{B?Ra)  in  Eq.  (A14)  is  the  approximation  for  the  most  singular  part  of  the  Green  s 
function  in  the  phi  integration  that  contains  the  self  patch.  For  the  remaining  less  singular 
part  of  the  Green’s  function,  can  be  approximated  by  the  mean  distance 

Ra  =  v^(At/2)V4  +  (p0')".  (A15) 

Combining  this  /?„  approximation  to  the  lesser  singular  part  of  the  Green’s  function  with 
in  Eq.  (A14),  one  obtains  Elq.  (47).  The  Rc  approximation  was  originally  derived 
for  R  so  that  the  phi  dependence  could  be  integrated  in  very  fine  increments,  but  after 
running  the  program  with  phi  integration  segments  much  finer  than  At,  it  was  determined 
that  phi-segments  approximately  equal  to  At  gave  the  greatest  accuracy.  It  was  also  seen 
that  the  accuracy  was  better  using  the  Rc  correction  than  the  correction  of  Reference  6. 
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ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re- 
searchf  development,  test,  and  technology  transition  in  support  of  Air 
Force  Command,  Control,  Communications  and  Intelligence  (C^l)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineertng  support  within 
areas  of  competence  is  provided  to  BSD  Program  Offices  (POs)  and  other 
BSD  elements  to  perform  effective  acquisition  of  C^I  systems,  bi  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


